High-order current correlation functions in Kondo systems. 
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We examine the statistics of current fluctuations in a junction with a quantum dot described by 
Kondo Hamiltonian. With the help of modified Keldish technique we calculate the third current 
cumulant. As a function of ratio v — bV/Tk the 3rd cumulant was obtained for three different 
regimes: Fermi liquid regime (y < 1), crossover interval (v > 1) and RG limit (v » 1). Unlike the 
case of noninteracting dot, 3rd cumulant shows strong non-linear voltage dependence. Only in the 
asymptotical limit v — > oo the linear dependence on V is recovered. 
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Motivation: The direct electron transport is an impor- 
tant tool to study small junctions and quantum dots. The 
investigation of the current correlation functions helps 
to get additional information about physical properties 
of such systems. During last years the measurements 
only of the second cumulant of fluctuating current (shot 
noise) have resulted in a remarkable theoretical and ex- 
perimental progress (see review article 0). Interest in 
the third and higher moments has occurred, first, be- 
cause its characteristics differ significantly from that of 
the second moment. In particular, as it was predicted 
by Levitov and Reznikov 0, in non- interacting systems 
the third moment is insensitive to the samples thermal 
noise, yet is more sensitive to the environment. Second, 
measurements of the higher moments may provide a new 
tool for studying conduction physics, complementary to 
the second moment. Experimentally it is quite nontrivial 
to extract cumulants higher then the second. Recently, 
however, the first experimental study of the third current 
cumulant in mesoscopic tunnel junctions was reported 
i. 

The theory of full counting statistic Q is a theoretical 
framework which is used to analyze statistics of charge 
transfer in experiment and to calculate the higher order 
cumulants However, when the interactions are in- 

cluded, which is the subject of the present paper, it is 
necessary to go beyond the theory of full counting statis- 
tics. Before starting such an analysis for Kondo type 
quantum dots we note that recently in a number of works 
third current cumulant in interacting mesoscopic systems 
was calculated 0, H| ■ In these works the 3rd current cu- 
mulat was defined as a time ordered product of three 
current operators on standard Keldish 2-time contour. 
As a result this correlation function is equal to a sum 
of partially time ordered products on the usual one line 
time contour. Thus we can apply Feynman perturbation 
theory to evaluate this correlation function, however, we 
are unable for interacting systems to compute non-time 
ordered 3rd cumulant. It is important to resolve this 
discrepancy. Here with the help of a modified Keldish 
technique which uses a multiple time contour (see Fig.l) 



we directly provide such a derivation of non-time or- 
dered 3rd cumulat for the Kondo problem. 

Hamiltonian. Third cumulant: We start with the 
Kondo Hamiltonian for quantum dot in the junction 
H = H L + H R + Hj where 

Hj = J2 ^*' c L(0)(^<W +Sw)c«v(0) (1) 

act' era' 

The first two terms correspond to non- 

1 

_ 2 .) 

(. 3 

4 ) 



FIG. 1: The four-axis time contour 

interacting electrons in the two leads -Hl(.r) = 
Y,k,a^L(R)kC S L{R)(jk c L (R) cr ^ where c a ^ k > £afc are 
the electron field operator and the electron energy of 
a lead. Index a = L. R indicates left (right) lead. We 
assume that the leads are dc-biased by applied voltage 
V. Here s is the one half spin matrix which acts on spin 
index of electron operators and S is the spin operator of 
the dot. The potential scattering is represented by the 
first term in the brackets. The bare coupling constants 
in Q can be obtained by Schrieffer- Wolff transformation 
from the parent Anderson Hamiltonian 

k,cr : a 

+ J2 ed l d a + Ud\d ] d\d l (2) 

cr 

and are related to the parameters of this Hamiltonian 0] 
as J aa < = v / r Q r Q '/(7r^|e|) and e = (U - \e\)e/U. Here 
U is the repulsive Hubbard coupling, e denotes bare level 
energy of the dot and v is the density of states in a lead. 
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We also introduce the widths Fl, R = 2itv\vl,r\ 2 which 
are expressed in terms of tunnelling matrix elements and 
density of electron states ^l). The current operator has 
a form 

1 = j[JLRY, C LM\^' + SSaa')c R ^(0) - H.C.] 

We define the 3rd momentum noise as a symmetric com- 
bination of non-time ordered correlation function of three 
currents. In a stationary situation this function can be 
written as 

s 3 (h - t 2 ,t 2 - 1 3 ) = I <i{t a )i{tp)i{t 1 )^z) 

P(a/3 7 ) 

where P(af3"/) is the permutation of 1,2,3. The Fourier- 
transformed value of £3(^1,^2) is a function of two 
energy variables. Below we consider the zero frequency 
limit £3(^1 = 0,ll>2 = 0) = £3 and introduce the third 
order cumulant S3, that is, the irreducible part of cor- 
relation function |3"|. This cumulant yields the equation 
S3 = S3 — 3/S2 + 2/ 3 . Here / is the averaged current and 
S2 stands for the pair current correlation function (shot 
noise). In the Kondo regime this function was recently 
calculated [l2|. 

Perturbation theory. RG: To apply perturbation the- 
ory to a product of non-time ordered Heisenberg oper- 
ators like S3 we need to order these operators on some 
multiple time contour. For two operators (shot noise) the 
standard 2-time Keldysh contour is sufficient. However, 
to arrange more then two current operators additional 
time axis must be included. The 3rd order correlation 
function S3 is describes by 4 lines contour (see Fig.l). 
9J. The action consists of the integration over this four 
lines path. However, as it is common for Keldysh tech- 
nique, we use only one infinite time path, though, with 
a four independent field operators, corresponding to the 
operators on each four branches ( Fig.l) of complete time 
contour. Thus we have: 

■ roo ^ 4 

S = di[£ cL.^iw - £ (4) 

J -°° kaa j = l 

Here & k , c aa ^ are the operators in 4(g) 4 Keldysh space, 

Gka denotes 4 ® 4 Green's function of noninteracting 
leads; a], = 1 if j=l,3 and a J z = —1 if j=2,4. The in- 
teraction now acquires Keldysh index j —>Hj- We 
use the important equalities that relate the new Green 
functions on the contour in Figl. to the ordinary Kcldish 
Green functions defined on two lines contour. Keeping 
explicitly only Keldish indices we can write them as: 

G 33 = G n , G 44 = G 22 (5) 
G»' = G 12 , j < j' 
= G 2 \ j>f 



On this four fold contour the correlation function is rep- 
resented by time ordered product of three current oper- 
ators 

S 3 (ti-t 2 ,t 2 -h) = \J2 <TP(t 1 )P(t 2 )I k (t 3 )> 

p(ijk) 

Here superscripts k are corresponding Keldysh in- 
dices of the Heisenberg operators (electrons and spin S) 
related to the first three paths in Fig.l (each index runs 
from 1 to 3). Summation over all permutations P of 
these indices is performed. Now we are ready to use 
non-equilibrium Keldysh technique and directly calculate 
S3. In the lowest (fourth) non- vanishing order of pertur- 
bation theory in the coupling constant J a . a ' the zero- 
frequency irreducible correlation function S3 acquires a 
form: S 3 = Sf ot + Sf + SJ*, with Sf ot = Y- 1 /(2h), 
Sf = W-x/H and S| x = 3Y 7 /(2h) where 

Y n = tlAlL J du(J h -f R ){f L (l-f R ) + 

f R {l-h)+nf L f R ) (6) 

Here Jl,b. are Fermi distribution functions for leads, 
and g aa i = J aa 'V stand for dimcnsionlcss coupling con- 
stants. To mark the origin of different contributions 
we intentionally separate S3 into three parts: the first 
one Sg * represents potential scattering, the second part 
Sg 6 defines mixed exchange and potential interactions, 
while the last term S% x stands for pure exchange scat- 
tering processes. For temperatures T >> Tk {Pk — 
D exp[— 1/(<?ll + g RR )], D Q are the Kondo temperature 
and the effective bandwidth, correspondingly) logarith- 
mic corrections which appear in perturbation theory can 
be disregarded. In this case the above formulaes describe 
S3 as a function of source-drain voltage (see inset in Fig. 2 
where S= S 3 ejPa\ and ob= Tte 2 g\ R /h~ is the conduc- 
tance in the Born approximation) for arbitrary values of 
the applied bias. It is interesting to note that unlike the 
case of noninteracting system a weak nonlinear voltage 
dependence of 3rd current cumulant at T » Pk persists 
in the lowest order of perturbation theory. 

The logarithmic divergences appear in the next (fifth) 
order expansion in couplings g aa i . We compute S3 to this 
order only for T = 0. In this case if voltages eV >> Pk 
the perturbation theory still can be used. We keep only 
those fifth order terms that consist of the maximal log- 
arithmic divergence. Only parts of S3 which include ex- 
change tunnelling are affected by these logarithmic cor- 
rections, while the pure potential part of the three cur- 
rents correlation function does not sufficiently changed. 
Thus after long, though, direct calculations we arrive at 

sr = ^_^V g A R[1 + 2(jgLL+gRR)bl Do ] (?) 
S? = ^f^9i R [l+4(9LL + 9 RR )hi^} (8) 
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When the voltage is decreasing the log terms are starting 
to increase so that for eV > Tk the expansions (|7I8() 
up to fifth order are not efficient. For this region 3rd 
momentum of the noise can be derived in the leading 
logarithmic approximation which consists of summation 
of the most diverging terms in each order in coupling 
constants g aa i . From a set of equations (|7I8() it becomes 
clear that the scaling behavior of S3 is similar to that of 
conductance 0. This becomes particulary clear when 
we rewrite the expression for the 3rd noise in a following 
form: 

S 3 (T = 0) = TiV^al + 3a B a (V) + \a 2 (V)\ (9) 

where a (V) = 3ne 2 [l + 2{g LL + g RR ) ]n(D /eV)]/(4K) 
is the 'spin exchange' part of conductance, the part 
which includes logarithmic term. On the level of 'poor 
man'fs scaling technique in the zero temperature limit 
the renormalization proceeds till the band width D be- 
comes equal to the applied bias eV. After that, the con- 
ductance can be calculated in the Born approximation 
with the renormalized exchange constant |lCj : o~(V) = 
37re 2 r L r jR ,/[2v / a(r L + r R )\n(eV/T K )} 2 . The potential 
scattering contribution to conductance ob is not changed 
under RG transformations and renders a small correction 
to S3. Thus, the final expression for the 3rd cumulant is 
given by equation Q where we should replace <jq(V) on 

<?{V). ' 

Fermi liquid regime: To study regime where T, eV < 
Tk we apply the mean field slave boson approximation 
(MFSB) |l3J. The current operator for Anderson model 
acquires a form 

T * = HE/[(^ c L(o)-««cL(o))^(o)-^.a] (10) 

In the slave-boson approach, the localized electron oper- 
ator d\ is represented by fib with b and fl being the 
standard boson and fcrmion operators. The total action 
on the time contour in Fig.l has the similar form as 
We should replace only Hj by interacting d-dependcnt 
part of the Anderson Hamiltonian @. The two new 
terms in the action, equivalent to the first one in equa- 
tion Q come from slave bosons and auxiliary fermions 
Green's functions. Also the requirement of a single occu- 
pancy b'b+Y^ a flfa = 1 must be included into the action 
with a Lagrange multiplier. In the MFSB approximation 
the Bose operators w,b are replaced by their expecta- 
tion value b. We also add to the action the source term 
S 7 = J dt"f 3 (t)J| (j=is running from 1 to 4) and inte- 
grate out the fermion operators in the leads. Here P s 
is the current operator for Anderson model taken on j 
branch of time contour in Fig.l. For a symmetric tun- 
nelling {vl = vji = v,Tl = T r = r ) the effective action 
can be written as S e ft(pf) — Sb + St where the first term 
is nonoperator bosonic part of the action, while the last 



one explicitly depends on source fields j k (i) and is given 
by 

Sfi-T) = zZj dt J dt 'tilGd*W]f° 
Gjfa) = G-^-niQ+GLvQ- + Q-G Ra Q+ 111) 

Here Gl,r represents 4 £g) 4 matrix of the electrons prop- 
agators in the leads, Gof a (u>) is the Fourier transform 
of zero order slave fermions Green's function and d z is 
diagonal 4 ® 4 matrix with elements u z 3 (@J. We also 
define T% — Tb 2 as an effective Kondo temperature. The 
G/<j(w) includes the Lagrange multiplier which shifts the 
localized level position e to e = e + A. Both free parame- 
ters Tfc and renormalized level e are self-consistently de- 
termined by two equations that define the extremum of 
S e f f relative to b and e when U — > 00 0] . 

We can trace out the slave fermions and get 
a closed form for generating functional: Z = 
expl—^j-TrlriGj 1 ^)]. The irreducible three currents 
correlation function or the 3rd cunmulant is obtained by 
taking the third variation of the logarithm of this func- 
tional on source fields. Performing the straightforward 
calculations we find S3 = S' 3 + S' 3 ' where 

ttp 3 r 

S' 3 = — j VM(/i-/fl) 3 (3-4^H) 

S * = IT / duj P 2 ^ L ~ Ml 1 - ^ L + f R - ^ 

Here the spectral density of interacting level irp(uj) = 
-ImG r d {u) = T k /[(uj - i) 2 + T 2 \ is introduced. If T -> 
then 3rd cumulant becomes particulary simple. Expres- 
sion for S'3 is now reduced to 

S 3 = -e 3 / dioT 2 (uj)(l - T(«)) (12) 
n Jo 

where T(u>)= nTkp(uj). Thus, we obtain 3rd noise S3 as 
a function of two dimensionless parameters eV/Tk and 
e/Ife. A simple approximate expression for S3 follows 
when these parameters are small S3 » 2e 6 V 3 j (3irT 2 ) . 
Such a voltage dependence is typical for a Fermi liquid. 

Discussion: The crossover regime (T, eV ~ Tk) actu- 
ally is extended to a several threshold values Tk- This 
region can also be studied within the present theory 
by using non-crossing approximation in non-equilibrium 
(NCA)[l4j together with 1/N expansion for generating 
functional Z . A satisfied estimation of the 3rd noise in 
this case can be achieved by calculating equation 112|) 
with a modified spectral density which is derived in NCA 
approximation: Tkp(uj) ~ * L / ojvca(w). We start numer- 
ical calculations by solving the mean field equations of 
MFSB theory. The MFSB method is limited to small 
voltages eV < Tk where is known it gives the correct 
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FIG. 2: The third cumulant S3 versus applied voltage eV/Tn- 
Dot curve represents MFSB calculations valid for small volt- 
age eV < Tk- Solid line with black squares is obtained in 
NCA approximation. Dash line is the result of direct com- 
putation of Eq.@ with ao — > o(V). The curves which were 
calculated for two Tk differing by a factor of 10 practically 
coincide. This reflects the universality of 3rd cumulant in 
considered regime, {inset): Normalized 3rd cumulant S_ as 
function of eV/T in the high temperature limit T >> Tk- 



qualitative behavior to this Fremi liquid regime. In addi- 
tion, at eV << Tk the 3rd cumulant in MFSB approxi- 
mation is matched with the perturbation theory deriva- 
tion at the unitary limit of Kondo Hamiltonian. For con- 
ductance and noise calculations 0| the whole scheme 
with MFSB works rather well. At eV < T K MFSB also 
gives adequate description of 3rd cumulant which is close 
to its exact value. The dot line in Fig. 2 displayed the re- 
sults of our MFSB calculations. 

In the crossover region the NCA computations begin 
with a finite V. The main block which we need to cal- 
culate in NCA is p(w). The 3rd current cumulant in- 
cludes the frequency integration of different powers of 
p(tjj). To estimate the accuracy of NCA we consider the 
known different physical values which are based on spec- 
tral density derivation. For example, for quantum dots, 
Cox 0] has shown that the calculated equilibrium sus- 
ceptibility agree with the exact Bethe ansatz results to 
within the 0.5% convergence accuracy of the NCA. For 
the current calculations Wingreen and Meir (see 0] ) use 
NCA equations in nonequilibrium and found at worse an 
overestimate of 15% on the linear response conductance. 
We expect the same (15%) accuracy for 3rd current cu- 
mulant. NCA result is given by solid curve with black 
squares in Fig. 2. Strong nonlincarity with a peak at the 
intermediate voltage eV m ~ 2.6Tk is determined by sig- 
nificant Kondo correlations. Indeed, in the absence of 
interacting S3 grows linearly with V. At eV ^ Tk we 
are in the strong Kondo regime where S3 grows consid- 
erably faster then simply linearly in V. The subsequent 



decreasing of 3rd noise for V > V m simply indicates on 
the weakening of the Kondo effect and manifests that 
Kondo correlations continues to determine the physics 
3rd noise. Even when we go outside the crossover region 
and enter scaling limit D n » eV » Tk (dash line on 
Fig. 2) where Kondo model is in a weak coupling regime 
with depressed Kondo effect, S3 shows a behavior which 
is far for being linear in V. Though, the Kondo effect is 
depressed, nevertheless, the logarithmic contributions to 
3rd cumulant which in this region are related to Kondo 
correlations are considerably larger than those derived 
with the help of bare perturbation theory in the Born 
approximation. Linear or weak nonlinear voltage depen- 
dence which would be indicate the complete depression of 
Kondo correlations appears outside the domain of Fig. 2 
for bigger T (as inset in Fig2) or V (T,V ~ D ). 

Conclusions: In conclusions, we have represented a 
general theory which allows to describe current fluctu- 
ations in quantum dot in the presence of interactions. In 
particular, we have considered a junction with quantum 
dot described by Kondo Hamiltonian. For this system 
we have calculated third cumulant of non-time ordered 
product of three current operators. Restricting ourselves 
to zero frequency we were able by our approach to cover 
three important regions: one is the in Fermi liquid limit, 
crossover regime, and weak coupling Kondo limit. Unlike 
the non-interaction system, in all cases S3 shows nonlin- 
ear voltage dependence. The linear dependence on V is 
restored only for large voltages v — > 00. In this work we 
assumed that external impedance is equal to zero, the 
condition, which may be violated in the real experiment. 
This impedance leads to a modification of the third cu- 
mulant |16| . 
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